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Abstract 

Consider a strong Markov process in continuous time, taking val- 
ues in some Polish state space. Recently, Douc, Fort and Guillin 
(2009) introduced verifiable conditions in terms of a supermartingale 
property implying an explicit control of modulated moments of hit- 
ting times. We show how this control can be translated into a control 
of polynomial moments of abstract regeneration times which are ob- 
tained by using the regeneration method of Nummelin, extended to 
the time-continuous context. 

As a consequence, if a p— th moment of the regeneration times 
exists, we obtain non asymptotic deviation bounds of the form 

f(X s )ds-n(f) 







Here, / is a bounded function and \x is the invariant measure of the 
process. We give several examples, including elliptic stochastic differ- 
ential equations and stochastic differential equations driven by a jump 
noise. 
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1 Introduction 



Let X be a positive Harris recurrent strong Markov process in continuous 
time, having invariant probability measure \x. From the Ergodic Theorem we 
know that for all x G 1R, f E L 1 (/i) and e > 



as t goes to infinity. The purpose of this paper is to establish the rate of con- 
vergence in (1.1), for bounded functions /. In the existing literature, mainly 
the case of exponential rate of convergence (exponential ergodicity) has been 
considered. But recently, there has been growing interest in studying other 
possible rates such as sub-geometric or polynomial rates. We follow this di- 
rection and study in this paper the case when the rate of convergence in (1.1) 
is polynomial. More precisely we use the so-called regeneration method and 
show that if a certain regeneration time admits a p— th moment, then we 
obtain non asymptotic deviation bounds of the form 



Here, / is a bounded function and [i is the invariant measure of the process. 
Such a bound is of major importance for many applications, for example non 
asymptotic problems for statistics of diffusions, concentration for particular 
approximations of granular media equations, and many other examples. 

Let us give some comments on the history of the problem and compare our 
result with known results on deviation inequalities for Markov processes. In 
the context of Markov chains, Clemengon (2001) and Bertail and Clemengon 
(2009) have obtained bounds in (1.1) which are exponential in time, using 
the regeneration method of Nummelin. They work under the conditions of 
geometric (exponential) ergodicity and stationarity, and within the space 
of bounded functions. Our work is close to this in spirit, since we use the 
regeneration method as well (however, we use it in a more complicated frame- 
work since we work in continuous time). Compared to their work we do not 
need to assume stationarity, our results hold for any starting point x or any 
starting measure provided it integrates the p— th moment of the regenera- 
tion time. Moreover, we weaken the assumption of exponential ergodicity 
to polynomial ergodicity. Still in the discrete framework of Markov chains 




(1.1) 




(1.2) 
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let us also mention the work by Adamczak (2009) who derives, using com- 
pletely different techniques, concentration inequalities for empirical processes 
of Markov chains, in the regime of exponential ergodicity. Finally, Chazottes, 
Collet, Kiilske and Redig (2007) obtain concentration inequalities for finite 
valued random fields on Z d via coupling both in the exponential and the 
sub-exponential regime. For their purposes, the finiteness of the phase space 
is crucial. 

All above mentioned results hold either in discrete time or in discrete 
state space, and this is not what we are interested in. In this paper we con- 
centrate on the framework of continuous time and general state space. For 
continuous time Markov processes there is a huge literature on the subject, 
and most of the results are based on functional inequalities and/or perturba- 
tion techniques which are used to obtain non-asymptotic bounds in (1.1). As 
a matter of fact, in contrary to our approach, most of these papers deal with 
the stationary case only or with the case when the initial law of the process is 
absolutely continuous with respect to the invariant measure, having a square 
integrable density. Wu (2000) uses the Lumer-Phillips theorem in order to 
derive non-asymptotic deviation bounds which are expressed in terms of a 
large deviations rate function. He works under the assumption that the ini- 
tial law of the process is absolutely continuous with respect to the invariant 
measure. Based on this, Cattiaux and Guillin (2008) use functional inequal- 
ities like the Poincare inequality in order to derive an exponential deviation 
bound; they work under the assumption of a spectral gap and with bounded 
functions. A small paragraph in Cattiaux and Guillin (2008) is devoted to 
the polynomial regime as well, under an assumption imposing polynomial 
decay of the a— mixing coefficient of the process, but the rate which is ob- 
tained is not optimal. In the same spirit, let us cite Guillin, Leonard, Wu 
and Yao (2009) who work in the space of Lipschitz functions under the as- 
sumption of a spectral gap. For bounded functions, they obtain a Hoeffding 
type inequality. Finally, Lezaud (2001) uses Kato's theory of perturbation 
of operators, still in the exponential regime. Let us also mention that in a 
completely different setting and having different applications in mind, Pal 
(2011) establishes concentration inequalities for diffusion laws on the path 
space C([0, oo)), using quadratic transportation cost inequalities. He stud- 
ies concentration around the median of the distribution, in the exponential 
regime, for Lipschitz functions on the path space with respect to the uniform 
norm. 

In contrast to most of the above mentioned papers, we do not assume 
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exponential ergodicity, nor the existence of a spectral gap nor stationarity. 
We do not need to assume that the process is \i— symmetric. The method 
we use is the so-called regeneration method. It appeals to the condition of 
integrability of regeneration times. Let us describe briefly what is the idea 
of regeneration times. In the easiest situation where the process X has a 
recurrent point x , we may introduce a sequence of stopping times R n , called 
regeneration times, such that 

1. For all n, R n < oo, R n +i = R n + Ri ° #R n , R n — > oo as n — > oo. (Here, 
9 denotes the shift operator.) 

2. For all n, X Rn = x . 

3. For all n, the process (X^ n+t ) t > is independent of J-j? n . 

In this case, paths of the process can be decomposed into i.i.d. excursions 
[Ri, Ri+i[, i = 1,2, ... , plus an initial segment [0, R±], and then limit theorems 
follow immediately from the strong law of large numbers. 

In general, recurrent points exist only in one-dimensional models. For 
one-dimensional recurrent diffusions it has been shown in Locherbach, 
Loukianova and Loukianov (2011) that, if for some p > 1 the p-th moment 
of the regeneration time exists, then (1.2) holds. 

For general multidimensional Harris recurrent processes, there is no direct 
way of defining regeneration times. However, there is a well-known method 
of introducing regeneration times artificially, which is known as method of 
"Nummelin splitting" in the case of Markov chains and which has been 
extended to the case of processes in continuous time by Locherbach and 
Loukianova (2008). This method consists of constructing a bigger process 
Z = (Z 1 , Z 2 , Z 3 ) taking values in E x [0, 1] x E, along a sequence of jump 
times = T Q < T x < . . . < T n < . . . , such that 

1. Z 1 is a copy of the original process X, and the T n are arrival times of 
a rate-1— Poisson process, independent of Z 1 . 

2. On each time interval [T n ,T n+ i[, Z 2 and Z 3 are constant. 

3. The sequence [Z\^) n is a copy of the resolvent chain X^ n+l (the process 
X observed after independent exponential times). 

4. The sequence (Z\ ) n is a copy of independent random variables, which 
are uniformly distributed on [0,1]. 
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The three co-ordinates and the sequence of jump times (T n ) n are constructed 
in a coupled way, inspired by the splitting technique of Nummelin (1978) and 
Athreya and Ney (1978) in discrete time. We recall the whole construction 
in Section 3. The main point of this construction is that there exist a mea- 
surable set C having pi(C) > (C will be a petite set in the Meyn-Tweedie 
terminology) and a parameter a G]0, 1] such that the successive visits of Z Tn 
to C x [0, a] x E induce regeneration times for the process Z. 

To resume, for any Harris recurrent Markov process X, the following 
holds true: the process X can be embedded as first co-ordinate into a new 
Markov process Z. This new process Z possesses regeneration times. These 
regeneration times are closely related to the hitting times of a certain petite 
set C, or in other words: the moments of regeneration times are closely 
related to hitting time moments. Once we have a p— th moment for the 
regeneration times, we obtain a control on the speed of convergence in the 
ergodic Theorem and (1.2) holds true. 

Note that different coupling techniques in spirit of the so-called Doeblin- 
or Dobrushin-coupling have been considered in the literature, for example in 
the case of diffusions by Veretennikov (1997) and (2004), and for Levy-noise 
driven solutions of SDE's by Kulik (2009). These couplings are more specific 
to the concrete models the authors are interested in - the coupling technique 
presented in this paper has the advantage of being completely general, as far 
as Harris recurrent processes are concerned. 

Once the coupling is constructed, it remains to establish sufficient condi- 
tions on the generator of the process ensuring that p— th moments for regen- 
eration times exist. These conditions are inspired by a recent work of Douc et 
al. (2009) on sub-geometric rates of convergence for strong Markov processes. 
In this work, the authors introduce a drift condition towards a closed petite 
set in the spirit of a condition of existence of a Lyapunov function. This 
condition provides an upper bound on the control of sub-geometric or poly- 
nomial moments of hitting times where the dependence on the starting point 
is precisely given. The drift condition also provides a verifiable condition 
ensuring positive Harris recurrence of the process. We recall these results in 
Section 2. Section 3 is devoted to give a self-contained description of the state 
of the art concerning the regeneration or Nummelin-splitting-method in the 
multidimensional case. Section 4 provides a link between the two approaches 
"Drift Condition" of Douc et al. (2009) and "Nummelin splitting". We 
show that the drift condition of Douc et al. (2009) provides an upper bound 
on the regeneration times introduced according to the method of Nummelin 
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splitting. More precisely, we show in Theorem 4.1 that certain polynomial 
moments up to a precise order are bounded - the bound on the order being 
determined by the Lyapunov condition. The dependence upon the starting 
point is controlled by the Lyapunov function as usual. So even though the 
moments of regeneration times can not be explicitly calculated, we get at 
least upper bounds in the rate of convergence in (1.1). As a main applica- 
tion of this result, in Section 5 we state and give the proof of the deviation 
inequality (1.2). Section 6 is devoted to some examples: multi-dimensional 
diffusions and SDE's driven by a jump noise that are treated in the spirit of 
a recent work of Kulik (2009). We close the paper with an appendix which 
recalls the Fuk-Nagaev inequality in the framework needed in Section 5. 

Acknowledgements. Eva Locherbach has been partially supported by an 
ANR project: Ce travail a beneficie d'une aide de l'Agence Nationale de la 
Recherche portant la reference ANR-08-BLAN-0220-01. 

2 Drift-condition, Harris-recurrence and mod- 
ulated moments 

Consider a probability space (Q,A,(P x )x)- Let X = (X t ) t > be a process 
defined on (Q,A, (P x ) x ) which is strong Markov, taking values in a locally 
compact Polish space (E,S), with cadlag paths. (P x ) x eE is a collection of 
probability measures on (Q,A) such that X Q = x P x — almost surely. We 
write (Pt)t for the transition semigroup of X. Moreover, we shall write (J 7 *)* 
for the filtration generated by the process. 

Throughout this paper, we impose the following condition on the transi- 
tion semigroup (Pt)t of X. 

Assumption 2.1. There exists a sigma-finite positive measure A on (E,£) 
such that for every t > 0, Pt(x, dy) = pt(x, y)A(dy), where (t, x, y) >->■ pt(x, y) 
is jointly measurable. 

We are seeking for conditions ensuring that the process X is recurrent in 
the sense of Harris. The most popular conditions for Harris-recurrence are 
drift conditions or more generally conditions in terms of a supermartingale 
property for a functional of the Markov process. We follow Douc et al. (2009) 
and impose a drift condition towards a closed petite set B which implies the 



6 



Harris recurrence of the process. Recall that a set B e £ is petite if there 
exists a probability measure a on B(M + ) and a measure z/ a on (E,£) such 
that 

/■oo 

/ P t (x,dy)a(dt)>l B (x)v a (dy). (2.1) 

Assumption 2.2. There exists a closed petite set 5, a continuous function 
V : E — > [1, oo[, an increasing differentiate concave positive function $ : 
[1, oo) —7- (0, oo) and a constant b < oo such that for any s > 0, x E E, 

E X (V(X S ))+E X {^J\oV{X u )dv^j < V(x)+bE x {^J\ B {X u )dv^j . (2.2) 

Remark 2.3. IfV & 'D(A) belongs to the domain of the extended generator 
A of the process X, then Theorem 3.4 of Douc et al. (2009) shows that 

AV{x) < -$ o V(x) + bl B (x) (2.3) 

implies the above Assumption 2.2. 

By Proposition 3.1 of Douc et al. (2009), we know that under Assumption 
2.2, the process X is positive recurrent in the sense of Harris. We write /x 
for its invariant probability measure. Hence, for any set A e £ such that 
(J>(A) > 0, we have limsup^^ 1a{XA = 1 almost surely. In particular the 
process is \i— irreducible. 

Under Assumption 2.2, Douc et al. (2009) give estimates on modulated 
moments of hitting times. Modulated moments are expressions of the type 

E x f r(s)f(X s )ds, 
Jo 

where r is a certain hitting time, r a rate function and / any positive measur- 
able function. Knowledge of the modulated moments permits to interpolate 
between the maximal rate of convergence (taking / = 1) and the maximal 
shape of functions / that can be taken in the ergodic theorem (taking r = 1). 
In the present paper we are interested in the maximal rate of convergence 
and hence we shall always take / = 1. 
For the function $ of (2.2) put 

f u ds 

H*(u) = ^ ^-y, u > 1, r*(s) = r(s) = $ o H^(s). (2.4) 
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We are interested in choices of the function $ that yield a polynomial rate 
function r. This is achieved by the choice $(i>) = cv a for < a < 1 giving 
rise to polynomial rate functions 

r(s) ~ Cs 1 ^*. 

We suppose from now on that Assumption 2.2 is satisfied with such a kind 
of function $(t> ) = cv a for < a < 1. The most important technical feature 
about the rate function that will be useful in the sequel is then the following 
sub-additivity property 

r(t + s) < c(r(t) + r(s)), (2.5) 
for t, s > and c a positive constant. We shall also use that 

r(t + s) < r(t)r(s), 

for all t, s > 0. 

We are interested in regeneration time moments. We will see in Section 3 
below that regeneration times are almost hitting times. Concerning hitting 
times, the following result is known in the literature. Fix 5 > and define 
for any closed set A e £ the delayed hitting time 

t a (5) := inf{t > 5 : X t G A}. 

Then Theorem 4.1 and Proposition 4.5, (ii) of Douc et al. (2009) imply the 
following two statements. Firstly, for the rate function r of (2.4) and for the 
petite set B of Assumption 2.2, 

r(s)ds < V(x) - 1 + J r{s)ds. (2.6) 

Second, for the rate function r of (2.4) and for any closed set A with fi(A) > 0, 
for any 5' > 0, 



E x [ Ai) r(s)ds < c(A, 5') \v(x) - 1 + / r(s)ds 

Jo L *v-U Jo 

Remark 2.4. Suppose that E = M and that the process X has continuous 
trajectories. Fix a recurrent point a G M. Then we can choose A = [a,oo[, 



(2.7) 
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if x < a, A —] — oo, a], if x > a in (2. 7) above. In this case, the successive 
visits 

Ri ■= T{a}(8), Rn+i ■= inf{t > R n + 5 : X t = a] 

of the point a are regeneration times of the process. Hence, (2.7) gives a 
control of regeneration time moments in the one- dimensional case. 

In the general multidimensional case, the times ta(S) do not define re- 
generation times any more. In this case, at least in general, regeneration 
times can only be introduced in an artificial manner, using the technique 
of Nummelin splitting in continuous time, as developed in Locherbach and 
Loukianova (2008). However, the estimates (2.6) and (2.7) can be translated 
into bounds on moments of these new extended regeneration times of the 
process. This is the main issue of this paper and will be treated in section 4 
below. 

In the next section we recall the technique of Nummelin splitting and then 
give the bounds on moments of the regeneration times. But before doing this 
we first recall some known facts about modulated moments of the resolvent 
chain from Douc et al. (2004). 

2.1 Modulated moments for the resolvent chain 

Observing the continuous time process after independent exponential times 
gives rise to the resolvent chain and allows to use known results in discrete 
time instead of working with the continuous time process. This trick is quite 
often used in the theory of processes in continuous time. 

Write U l (x,dy) := / °° e~ t P t (x ) dy)dt for the resolvent kernel associated 
to the process. Introduce a sequence (<7„)n>i of i.i.d. exp(l)-waiting times, 
independent of the process X itself. Let To = 0, T n = o\ + . . . + o~ n and 
X n = Xt„- Then the chain X = (X n ) n is recurrent in the sense of Harris, 
having the same invariant measure \i as the continuous time process, and its 
one-step transition kernel is given by U x {x, dy). 

Since X is Harris, it can be shown (Revuz (1984), see also Proposition 
6.7 of Hopfner and Locherbach (2003)), that the resolvent satisfies 

U\x,dy) > al c (x)u(dy), (2.8) 

where < a < 1, n{C) > and v a probability measure equivalent to 
//(• fl C). The set C is in general not the petite set of Assumption 2.2. It 
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can be chosen to be compact. In particular, (2.8) implies that the resolvent 
chain is aperiodic. 

It is interesting to note that the drift condition (2.2) on the process in 
continuous time implies a similar drift condition on the resolvent chain. More 
precisely, Theorem 4.9 of Douc et al. (2009), item (ii), implies that under 
Assumption 2.2 the resolvent chain satisfies a drift condition as well, with a 
different petite set and different functions $ and V, but giving rise to the 
same rate function r since $(£(1 + $'(1))) ~ $(t) for t — > oo. Moreover, 

||t7-y(i + $'(i))|| 00 <oo. 

Now for any measurable set A with fi(A) > 0, write ta '■= inf{n > 1 : X n G 
A}. Then, by Douc et al. (2004), proof of Theorem 2.8, second formula, 



TA-l 



k=0 



< c^Afrix) + c 2 {A) < Cl V(x) + c 2 , 



(2.9) 



since V(x) < CiV(x) + c 2 . 

After these preliminaries on resolvent chains we now turn to the descrip- 
tion of the regeneration method in the case of a general state space. 



3 Nummelin splitting and regeneration times 

Regeneration times can be introduced for any Harris recurrent strong Markov 
process under the Assumption 2.1 - without any further assumption. We 
make once more use of the resolvent chain. Recall the definition of the 
resolvent kernel U 1 and the lower bound (2.8) which holds under the only 
assumption of Harris recurrence: 

U l (x,dy) > al c (x)v(dy), 

where C is a fixed compact petite set with /i(C) > 0. Note that since [x(C) > 
0, (2.7) and (2.9) hold for the hitting time of this set C. 

Remark 3.1. Fort and Roberts (2005) and Douc et al. (2009) impose quite 
systematically the condition of irreducibility of some skeleton chain, see e.g. The- 
orem 3.2 and Theorem 3.3 of Douc et al. (2009). This implies the existence 
of some m such that P m satisfies 

P m (x,dy) > al c {x)v{dy). 
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This condition is obviously stronger than (2.8) and implies that the process 
is not only positive Harris recurrent but also ergodic, i.e. for all x e E, 

\\P t (x, .) - h\\tv -> 0. 

We do not impose this additional condition. 

We now show how to construct regeneration times in continuous time 
by using the technique of Nummelin splitting which has been introduced 
for Harris recurrent Markov chains in discrete time by Nummelin (1978) and 
Athreya and Ney (1978). The idea is to define on an extension of the original 
space (Q,,A,(P X )) a Markov process Z = (Z t ) t > = (Z} , Z?, Zf) t > , taking 
values in E x [0, 1] x E such that the times T n are jump times of the pro- 
cess and such that ((Zl) t , (T n ) n ) has the same distribution as ((X t )t, (T n ) n ). 
We recall the details of this construction from Locherbach and Loukianova 
(2008). 

First of all, define the split kernel Q{{x, u), dy). This is a transition kernel 
Q((x, u), dy) from E x [0, 1] to E defined by 

(vidy) if (x,u) eC x [0,a] 

{U\x, dy) - av(dy)) if (x, u) e Cx]a, 1] (3.1) 
U^x.dy) iixiC. 

Remark 3.2. This kernel is called split kernel since duQ((x,u), dy) = 
C/ 1 (a;, dy). Thus Q is a splitting of the resolvent kernel by means of the addi- 
tional "colour" u. 

Write m 1 (x, x') := J* °° e _ *p t (x, x')dt. We now show how to construct the 
process Z recursively over time intervals [T n ,T n+1 [,n > 0. We start with 
some initial condition Z\ = X = x, Zq = u G [0, 1], Z$ — x' e E. Then 
inductively in n > 0, on Z Tn = (x, u, x') : 

1. Choose a new jump time a n+ i according to 

, p t (x,x') , m 
e i / a dt on R+, 
u l (x, x') 

where we define 0/0 := a/oo := 1, for any a > 0, and put T n+i := 
T n + a n+ i. 

2. On {<7 n+ i = t}, put Z\ n+S := u, Z\ n+S := x' for all < s < t. 
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3. For every s < t, choose 



! p a (x,y)p t -s(y,x') w, x 
^ P. (*,*') A( ^ 



Choose Zr n+S := x for some fixed point x & E on {p t (x,x') = 0}. 
Moreover, given Zj, +s = y, on s + u < t, choose 

7 i Pu(y,y')Pt- s -u(y',x') 
Z ™ u pT^) m) - 

Again, on {p t - a (y, = 0}, choose Zl +S+U = x . 



4. At the jump time T n+1 , choose Z\ n l := Z\ n = x' . Choose Z^ n 1 inde- 
pendently of Z s , s < T n+ i, according to the uniform law U. Finally, on 
{Zl +i = «'}, choose Z* n+i ~ Q((x', «'), tfo"). 

Note that by construction, given the initial value of Z at time T n , the 
evolution of the process Z 1 during [T n ,T n+1 [ does not depend on the chosen 
value of Z% . 

We will write P n for the measure related to X, under which X starts 
from the initial measure ir(dx), and IP^ for the measure related to Z, under 
which Z starts from the initial measure n(dx) ®U(du)®Q((x, u), dy). Hence, 
IP Xo denotes the measure related to Z under which Z starts from the initial 
measure 5 X(j (dx) <8> U{du) <8> Q((x,u),dy). In the same spirit we denote E n 
the expectation with respect to P w and JE n the expectation with respect to 
lP n . Moreover, we shall write JF for the filtration generated by Z, G for the 
filtration generated by the first two co-ordinates Z x and Z 2 of the process, 
and W x for the sub-filtration generated by X interpreted as first co-ordinate 
of Z. 

The new process Z is a Markov process with respect to its filtration IF. 
For a proof of this result, the interested reader is referred to Theorem 2.7 of 
Locherbach and Loukianova (2008). In general, Z will no longer be strong 
Markov. But for any n > 0, by construction, the strong Markov property 
holds with respect to T n . Thus for any f,g:Ex [0,1] x E — >■ 1R measurable 
and bounded, for any s > fixed, for any initial measure n on (E,E), 

EMZ Tn )f(Z Tn+s )) = EMZ Tn )lE ZTn (f(Z s ))). 

Finally, an important point is that by construction, 

C{{Zl) t \P x ) = C{{X t ) t \P x ) 
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for any x e E, thus the first co-ordinate of the process Z is indeed a copy 
of the original Markov process X, when disregarding the additional colours 

[y. 1 . ■/■'). 

However, adding the colours (Z 2 ,Z 3 ) allows to introduce regeneration 
times for the process Z (not for X itself). More precisely, write 

A := C x [0, a] x E 

and put 

S := 0, R := 0, S n+1 := inf{T m > R n : Z Tm E A}, 

R n+l := inf{T m : T m > S n+1 }. (3.2) 

The sequence of IF— stopping times R n generalises the notion of life-cycle 
decomposition in the following sense. 

Proposition 3.3. [Proposition 2.6 and 2.13 of Locherbach and Loukianova 
(2008)] 

a) Under IP X , the sequence of jump times (T n ) n is independent of the first 
co-ordinate process {Z]) t and (T n — T„_i) n are i.i.d. exp{V) — variables. 

b) At regeneration times, we start from a fixed initial distribution which does 
not depend on the past: Zr u ~ u(dx)U(du)Q((x,u),dx') for all n > 1. 

c) At regeneration times, we start afresh and have independence after a wait- 
ing time: Z Rn+ . is independent of Fs n - for all n > 1. 

d) The sequence of (ZR n )„>i is i.i.d. 

Since the original process X - under Assumption 2.2 - is Harris with 
invariant measure fi, the new process Z will be Harris, too. We shall write 
n for its invariant probability measure, n can be written in terms of an 
occupation time formula which is a consequence of Chacon-Ornstein's ratio 
limit theorem. In order to state this theorem, let us recall that an additive 
functional of the process Z is a M + — valued, IF— adapted process A = (A t ) t > 
such that 

1. Almost surely, the process is non- decreasing, right-continuous, having 
A = 0. 

2. For any s,t > 0, A s+t — A t + A s o9 t almost surely. Here, 9 denotes the 
shift operator. 
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The additive functional is called integrable if E u (Ai) < +00. Examples for 
integrable additive functionals are A t = J Q f(Z s )ds, where / is a positive 
measurable function, integrable with respect to the invariant measure II. 

Proposition 3.4 (Chacon-Ornstein's ratio limit theorem). Let A t , B t be any 

positive additive functionals of Z such that En(Bi) > 0. Then 

A t JEufAA m 

— > , . JP X — almost surely, as t — > 00, 

B t JE n {Bi) 

for any x e E. Moreover, Z is recurrent in the sense of Harris and its unique 
invariant probability measure II is given by 

11(f) = £E W f(Z s )ds, (3.3) 

JRi 

where t = E(R 2 - Ri)~ l > 0. 

Proof. The proof follows easily from the regeneration property with respect 
to the regeneration times R n . □ 

The invariant measure fi of the original process X is the projection onto 
the first co-ordinate of II. From this we deduce that the invariant probability 
measure \x of the original process X must be given by 

rRi 

n(f)=eiE n f(X s )ds, (3.4) 

JRi 

where we recall that I = E(R 2 — Ri)' 1 > 0. In the above formula we 
interpret X as first co-ordinate of Z, under lP n 1 . R 2 — R\ is the length 
of one regeneration period. Under assumption (2.2), the process is positive 
recurrent and hence the expected length t of one regeneration period is finite. 

We now turn to the main issue of this article which is the control of the 
speed of convergence in the ergodic theorem. As a consequence of the above 
considerations, we can write 



Px 



- f f(X s )ds-fi(f) 
1 Jo 



> 5 



1 rt rn.i 

- \ f(zl)ds-eE„ / f(z])d 

1 Jo JRi 



R 2 

" > 5 ) , (3.5) 



1 Actually, we should write f^ 2 f(Zl)ds - but if not otherwise indicated, this iden- 
tification will always be implicitly assumed. 
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where we recall that JP X denotes the measure related to Z under which Z ~ 
5 X <g> U(du) <S> Q((x,u),dy). The more moments of the regeneration period 
R2 — R1 exist, the more the process is recurrent and the more the convergence 
in (3.5) is fast. 

We first give estimates on the polynomial moments 




depending on the starting point x. Integrating this against v(dx) gives then 
a control on the corresponding moment of the regeneration period. This 
integration does not pose any problems because the support of the measure 
v is the compact set C. Since our regeneration times are built based on the 
resolvent chain, the main technical ingredient that allows such a control will 
be the estimate (2.9) rather than (2.7). 



4 Polynomial moments of regeneration times 

The aim of this section is to show that the results of Douc et al. (2009) 
can be translated immediately into a control of moments of regeneration 
times. This yields somehow a link between the two different approaches 
"Drift conditions" versus "Nummelin". Recall the definition of r(s) = r$(s) 
in (2.4). 

Theorem 4.1. Grant assumptions 2.1 and 2.2 with a function $(t>) = cv a , 
where < a < 1. Then there exist constants c\ and c 2 , such that 

rRi 

E x I r{s)ds < c 1 V(x) + c 2 . 
Jo 

Remark 4.2. For $(v) = cv a , it can be easily shown that there exists a 
constant c such that r(s) = r$(s) > cs~. Hence the above theorem implies 
the control of polynomial moments of the regeneration time, i. e. 

E x Rj^ <c 1 V(x) + c 2 . (4.1) 
Proof. Recall the definition of the regeneration times in (3.2). Let 
5i := inf{T n : Z\ n G C}, S n+l := M{T k > S n : Z l Tk G C}. 
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Obviously, R\ > Si. 



1. In the following, c will denote a constant that might change from line 
to line. We first show how to control 

E, 



f 

/ r(s)ds. 
Jo 

In a first step we show that 

rS\ poo poo 

E x r(s)ds = E x e~ti lc{zl)ds r(t)dt = E x e~ f ° lc{x ^ ds r(t)dt. 
Jo Jo Jo 

(4.2) 

This can be seen as follows. First, in order to obtain the law of Si, we 
evaluate for any a > 0, 

P x (Si>a) = Y, F -CSi = T n ,T n >a) 

n>l 

= J2 F *( z k e C c , . . . , z\ n _ x e C c , z l Tn eC,T n > a) 

n>l 
n>l 

= E x ( £(1 - l c (X Tl )) • • • (1 - 1 c (X Tb _ J)/(X Tb , T B ) J , 

\n>l / 

where = l t>a l c (x). 

Now, we make use of the following very useful formula which is taken 
from Hopfner and Locherbach (2003), (5.29), page 59. 

E x ( £(1 - l c (X Tl )) • • • (1 - 1 C (X T „_J)/(X T „, T n ) J 

\n>l / 

= £ x QH l c (X t )e"-tf M*.)* d ^ . 
Hence we obtain 

PsOSi > a) = E x QH l c (X t )e-J"o lo ^) <fa d^ = E x ( e -.fo M*.)*) 
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Writing finally that 



po\ poo poo 

E x / r(s)ds = E X l s<§i r(s)ds = / r(s)F x (5i > s)ds, 
Jo Jo Jo 

we get (4.2). No we apply once more formula (5.29) of Hopfner and Locherbach 
(2003) and obtain 

roc / 00 \ 

E x J e fo ic(x s)d s rm = Ex I ^(1 - \ C (X X )) ■■■{!- l c (X B _ 1 ))r(T B ) j , 

(4.3) 

where we recall that X n = X Tn is the process observed at the n— th jump 
time of an independent rate one Poisson process. The expression at the right 
hand side of (4.3) is almost a modulated moment for the resolvent chain, 
except that we have to replace r(T n ) by r(n). This is not difficult since for n 
large we can use the law of large numbers. Since r is increasing we can write 

E x ((1 - 1 C (X0) • • • (1 - l c (X n _ 1 ))r(T„)) (4.4) 
< E x ((1 - • • • (1 - l c (X n _0)r(2n)) 

+E X ((1 - l^XO) • • • (1 - lc(X„-i))lT„> 2 „r(T n )) .(4.5) 

Let us start with the control of the first term in this decomposition. Recall 
that fc = inf{n > 1 : X n G C}. Now, using that r(2n) < cr(n), which 
follows from r(t + s) < c(r(t) + r(s)) by (2.5), 



E x ]T(1 - 1 C (X0) • • • (1 - l c (X n _ 1 ))r(2n) = ^ £ r(2r 

\n=l / \n=l 

^ r(n) J < c£?J ^ r(ra) J + cE x r\f c ). (4.6) 

Let = Xlj=o r (i)- Since r is polynomial, lim^oo r(k)/R(k) = 0. Hence 
there exists a constant c such that for all k > 1, r(fc) < + c. As a 

consequence, 

/rc-l 

^r(f c ) < c + 1E X r ( n ) 

V n=0 

Using (2.9), we can thus conclude that 

(oo 
- M*0) • • • (1 - l c (X n _i))r(2n) ) < d^x) + c 2 . 
n=l 
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Now we turn to the second expression in (4.4) above: For any 1 < p, q such 
that i + i = 1, 

E x ((1 - lcpTO) • • • (1 - lc(^n-l))lT n >2r»r(T n )) 

< [iE^(T n )] 1/p -[iP,(T n >2n)] 1/( ' 

< [E x r^T n )\^ ■ e~ Cn 



(4.7) 



for some suitable constant C. But r p (-) is polynomial and T n the sum of n 
independent exp(l) variables, hence sup x IE x r p (T n ) < P(n), where P(.) is a 
polynomial in n. As a consequence, 

^ sup iE, ((1 - ldXj) • • • (1 - l c (X„_i))l Tn>2n r(T„)) = C 2 < oo. 

n>l X 

Putting together (4.2), (4.3)-(4.7), we thus get that 



E x ! r(s)ds < CiV(x) + c 2 . 



(4.8) 



This will be the main contribution to the control of 1E X J Q Rl r(s)ds. In the 
sequel, we shall also use that (4.8) implies in particular 



/•Si 

sup 1E X / r(s)ds < +oo, 
xec Jo 



(4.9) 



since C is compact. 

2. Recall the definition of Si in (3.2). We now show how to use the control 
of Si in order to obtain a control of Si. We have, since r(t + s) < r(s)r(t), 



IE, 



Si 

, / r(s)ds 
o 



jEr. 



Si 



fSl I rSn+l—Sn 

E x / r(s)ds + ^E X \ / r(S n + s)dsl §n<Sl 

J n >l \J0 

f 

/ r(s)ds 
Jo 

I rS n +i—S„ 

+ / r(s)ds r(Sn)l,s n<S i )• (4-10) 



n>l 



Sn-i 



S )ds + L r ( S ) dsl S n <Si 



n>l 





l"S„+l—S n 




( 


/ r(s)ds 


r (^) 1 5„<5 1 ^ 




Jo 
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The first term in this expression can be controlled using (4.8). We study the 
second term in the above expression 



®* \r(S n )ls n<Sl j Sn+1 S r{s)ds 



o 



We know that lP x (S n < Si) — (1 — a) n (see for example the proof of Propo- 
sition 2.16 in Locherbach and Loukianova (2008)). A first idea would be to 
use Markov's property with respect to S n : 



/ _ pSn+i — S„ \ I /-Si 

®* [r(S n )ls n<Sl j o r(s)ds \ = E x I r(S n )l §n<Sl lE z§n ^ r(s)ds 



But unfortunately it is not true that 

"§! /-Si 



JEz- I r(s)ds < sup E x I r(s)ds, 
Sn Jo xec Jo 

we only have that on {S n < S±}, 

rSi pSi 

JEzg / r(s)ds < sup lE (x ^ z) / r(s)ds, 

n JO x£C,u>a,z€E Jo 

and this can not be directly controlled using (4.8). 

Hence, we must be more careful. We use that r(S n )l^ n<Si y is measurable 
with respect to Qg where we recall that (Gt)t is the filtration generated by 
the first two co-ordinates Z 1 and Z 2 of Z. Hence we will condition on Qa ■ 
Note that by construction of Z, this means that we condition on the whole 
history of the whole process, i.e. the three co-ordinates, up to the last jump 
time sup{Tfe : Tk < S n } strictly before S n , and on the history of Z 1 and Z 2 
up to time S n . In other words, conditioning on , we know Z\ and Z\ , 

but Z% has still to be chosen. Moreover, on {S n < SA, Z\ > a, and hence 
the second line in the definition of the kernel Q((x, u), dx') of (3.1) has to be 
applied. 

Write v{x) for the density of v{dx) with respect to the dominating mea- 
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sure A(dx) of assumption 2.1. Then, 

'o 



(Sn)l§ n<Sl j o r ( S ) ds 

'' " \Zl n ,x')-au(x')j A(dx') 



/•oi 

< y^E x (r(S n )ls n<Sl J E V>{Z\ n ,x')\{dx') 

E (z , >z| ^ / \{s)ds] .(4.11) 



But for any u, 



E 



u (x,x')A(dx')]E^ UjX ,) / r(s)d 



Si 



= du Q((x,u),dx')E {x ^ x/) / r(s)ds, (4.12) 
Jo Je Jo 

since JE^ X)UtX ^ J^ 1 r(s)ds does not depend on u. Moreover, 

/■l /■ /-Si />Sl 

du Q{{x,u),dx')lE( x ^ x >) / r(s)ds = E x I r(s)ds. 
Jo Je Jo Jo 



Hence, since Z\ e C, 



^ I r(5 ri )l^ <5i / r(s)ds 



(4.13) 



< -^—E x \ r{S n )ls n<Sl ^sup B x r(s)ds 

< Y^E* (r(S n )ls n<Sl ) ■ (4.14) 
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Hence we must control E x (lg n<Si r(S n )). We write S n = S 1 + (S n - S 1 ) 
and use once more the sub-multiplicativity of r. We obtain 

E* (<Sn)ls n<Sl ) < E x (riS^sASn - Si)l §n<Sl ) . (4.15) 

Here, we have cut S n — Si + (S n — Si) into two pieces in order to get a last 
term which does not depend on the starting point. The same arguments as 
above in (4.11) and (4.12) yield, when conditioning on Q§ , the following. 

E* (r(S n )ls n<Sl ) < E x (r(Si)ls 1<Sl r(S n - Si)ls n < Sl ) 
<E x (r(Si)ls 1<Sl ^Ju\Zl,x')A^ 

(4.16) 

Concerning the last term in the above expression, we use that r(t) < J* * r(s)ds+ 
c for some constant c and obtain 

E x (r{Si)ls 1<Sl ) < c + E x (^J^ r(s)ds 

< c + ciV(x) + c 2 = ciV(x) + c 2 , (4.17) 

using (4.8). 



Concerning the first term in (4.16), for p, q > 1 such that ^ + ^ = 1, we 
obtain 

[E u r"(S n . L )) 1 ''" 

yeC ^ ' y&C 



supE y (r(S n -i)ls n 1<Sl ) < sup(E y r p (S n -i)) * E y (S n „i < Si) 1 ^ 
y ec ^ ' yec ^ / 

< (1 -a^-^f sup E y rP(S n _i))) %.18) 
\yec J 

We have to control this last expression. We claim the following: There exists 
a k > and a constant c such that for p > 1 sufficiently small, 

supE y r p (S n -i))) <cn K . (4.19) 
v&c J 

21 



Once (4.19) is proven, we obtain, using (4.10), (4.13), (4.16), (4.17), (4.18) 
and (4.19) the following: 

Si 



r(s)ds 



< ( Cl V(x) + c 2 ) + 



{c 1 V(x) + ~c 2 )Y J ^-a) (n - l),q n K 



n>l 



= c^Vix) + c 2 . (4.20) 

It remains to show (4.19): By our assumptions, r is polynomial and 
r(x) ~ Cx 1 ^ as x — > oo, hence r p (x) < cx Kp , where n = a/(l — a). We now 
fix the choice of p and q in (4.18). We choose 

,1 1 r . 1 — a 1 r 
pe]-l + -[ = ] , - . 

k k a a 

Then Kp > 1, and we use Jensen's inequality to obtain 

r p (S n ^) < cSZi <(n- l)""" 1 fe p + . . . + (5 n _i - 5„_ 2 ) Kp 



(4.21) 

Now since p < 1 + ^ = ^, we have t Kp < c f*r(s)ds for some constant c. 
Then for any of the above terms (k > 2), by (4.9), 



sup lE y (Sk — Sk-i) Kp < csupiE y / r(s)ds < oo. 



As a consequence, coming back to (4.21), 



Si 



sup Ey^Sn-x) < c(n - l) pK sup ^ / r(s)ds = c(n - l) pK , 
yec y£C Jo 

and this yields (4.19). 

3. Finally we proceed to the control of R\. Clearly, 



E 



rRi 

Jo 



Si 



r(s)ds < E x / r(s)ds + E x 



Ri-Si 



r(s)ds 



We have to control the last term above. We condition on Qs 1 , notice that 
Zsi — a an d use s ^ e P 1- °f the construction of Z, hence 



r(Si) 



Ri-Si 



= JE X 



r(s)ds 
u(x')A(dx') 



E 



' P i Z k X 'U t fr{s)ds 
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But by (2.8), v{x') < ^(Z^x'), since Z\ n e C, thus 



E, 



r(5i) / 

Jo 



r(s)ds 



< -IE,. 



a 



r(^) ( f A(dx') e-tptiZ^x^dt [ r(s)ds 
\Je Jo Jo 



= —lEr, 



a 
1 

a 



r(5x) ^ e-'dtlf Pt (Z l Si ,x')K{dx')\[f r(s)d 

Jo \J E J Jo 



-E^ 



r(S 1 ) / e-*dt / r(s)ds 



-^(r^)), 



since J °°e * J* * r(s)dsdt < 00. Finally, r(t) < J^r(s)ds + c gives 



))<JE X 

Jo 



r(s)ds + c, 



which is controlled due to (4.20). This concludes the proof. 



□ 



Remark 4.3. The fact that the rate function is polynomial was crucial at 
two points in the above proof: in equations (4-7) and (4-19). The general sub- 
geometrical case could probably be handled by paying in particular attention 
to the constants that arrive in expressions like E x r p (T n ) < [E x r p (Ti)] n . 



5 Polynomial deviation inequality 

We impose Assumption 2.2 with a function $(t>) = cv a , where < a < 1. 
As a consequence, we obtain a control for polynomial moments E X R\ of the 
regeneration time for all p < 1/(1 — a), due to (4.1). Since V is continuous 
and since the measure v of (2.8) which is used in order to construct the 
regeneration periods is of compact support, also E V R\ is finite for all p < 
1/(1 -a). 

In order to derive the deviation inequality we first state a deviation in- 
equality for the counting process associated to the life cycle decomposition 

00 

N t = sup{n :R n <t} = J2 l {R n <t}, N = 0. 

n=l 
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We have almost surely, as t — > oo, N t /t — > JE u Ni = £, where we recall that 

£ = {E^Y 1 = {E{R 2 - RJ)- 1 , 

see Proposition 3.4 and equation (3.3). 

The deviation inequality for the counting process associated to the life 
cycle decomposition is the following. 

Theorem 5.1. Grant Assumptions 2.1 and 2.2 with $(t>) = cv a , < a < 1. 
Let x G E be any starting point and < e < 1. Then for any 1 < p < 
1/(1 — a) there exists a positive constant C(£,p,u) such that the following 
inequality holds: 

** (if - il > i£ ) ~ C «> P > V)V{X) gPAg Vi) ^T" ^ 



Here C(£,p, u) is given by 
C(t,p,u) = 



C{p) 
C(p) 



(1 + (l/ey- 1 ) + (m p + a 2 ^)] (IP' 1 V £) ifp > 2 
(l + (l/£y- l )+m p tCv\ i/pe]l,2[ 



where C(p) is a constant depending only on p, C p is the constant of the 
Burkholder- Davis- Gundy inequality, m p = E(R 2 —Ri — j) p and a 2 = Var(R 2 - 
i?x), in the case p > 2. 

Proof. The proof is basically the same as in Locherbach, Loukianova and 
Loukianov (2011), proof of Theorem 3.1. Put in contrary to there we use 
the Fuk-Nagaev inequality given in the appendix (Theorem 7.1) in the case 
p > 2. We decompose: 

P x (\N t /t -£\> is) < P x (N t /t > £(1 + e)) + P x (N t /t < £(1 - e)) . (5.2) 

Put for k>l,fjk — —l(Rk+i — Rk — 1/0- For the fi rst term °f (5-2), we 
have 

/ [U(l+e)\ 

P x (N t /t > £(l + e)) < P x (Rx - l/£ < -te/2) + P x Yl Vk> te/2 

y k=i 

(5.3) 
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In an analogous way, 

, l \ / [U(l-e)]-l 

JP X (N t /t < i{\ -e))<P x \R x -j> te/2 J + iPj £ jfc < -te/2 

(5.4) 

The random variables f]k,k > 1, are identically distributed centred random 
variables such that IE x \fj k \ p < oo. Moreover, they are two-dependent. Indeed, 
fjk is not independent of J-'r,,, but only independent of J r R k _ 1 - This is due to 
step 1 of the construction of Z, where the waiting time for the new jump is 
chosen depending on the actual value of Z at time R k - 

If p > 2, we can apply Theorem 7.1. Let M = and M n = Ylk=i 'fa- 
Denote M* = sup fc<n \M k \. As a consequence of (5.3) and (5.4) we can write 

P x (\N t /t -£\> is) < P x (IR, - l/£\ > te/2) + P x (M{ m+e)] > te/2) . 

(5.5) 

We use Theorem 7.1 with n = [t£(l + e)] and A = te/8 and obtain 
P x (\N t /t-e\>£e) 

< 2P " 1 ^ 1 p ~ 1/ ^ |P " 1 + C(p)[m p + a 2 ^] (gr 1 V i) e -2(P-i) r (P-i) 

< (2^ 1 ^|i? 1 - 1/^r 1 + C(p)[m p + a 2 ^" 1 )] (F- 1 V €)) -J-^r^" 1 ), 

since e < 1, where m p = lE x \fji\ p , o 2 = Var(fji). Finally we use that 

E X \R, - l/ir 1 < C{p)[lE x R\- 1 + (l/iy-'i 
and that for some constants c and d, 

JE x {R p ~ l ) < 1 + E X R\ < cV(x) + d 
to conclude that, since V(.) > 1, 

P x (\N t /t -£\> is) < C(p) V(x)((l + (l/£) p ~ 1 )+ 

K + ^)] (f-v^))^-^). 

This finishes the proof in the case p >2. 
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In the case 1 < p < 2, we apply the Burkholder-Davis-Gundy inequality. 
In order to produce independent random variables, we define 

(i) _ f fjk if A; odd 1 (2) _ f fj k if fc even 1 
^ _ \ elseif /' Vk ~ \ elseif J ' 

Let Mo 1 = and = ELi 7 ?/^- In the same W, = and M l = 

Efc=l ^ ■ 

We also introduce the following two sub-filtrations, associated to the sum 
of odd and the sum of even terms. Let 

:= a{rff : k < n,k odd } = a{M%\k < n}, 

and 

:= a{ri ( ^ :k<n,k even } = a{M^ 2 \k < n}. 

Then (M^) n and (M 2 )„ are discrete A& —martingales (An''— martingales, 
respectively). Both martingales are L p martingales such that [MW]„ = 

ELifa* ) 2 . for i = !.2- Denote ( M(l) )n = sup fc < n |M«|, i = 1,2. In this 
case, as a consequence of (5.3) and (5.4) we write 

P x (\N t /t -e\> is) < P x (|i?i - > te/2) (5.7) 
+P X ((M<% (1+e)] > te/4) + F x ((M( 2 ))f« (1+£)] > te/4) . 

We use the Burkholder-Davis-Gundy inequality to bound the last terms 
in (5.7): For all p > 1 there exists a constant C p depending only p such that 

ll(M«);|| p < C P \\[M%%, hence ^((JlK^ < C>E X feLrf,^'' " 



Notice that by definition, the term Efe=i(% ) 2 contains [^^] terms 
whereas Efc=i( r ?i 2 ' ) ) 2 contains [ra/2] terms. Since 1 < p < 2, the sub- 
additivity of the function x >->■ x P//2 implies 




fc=i 

The same kind of bound holds also for the even terms. 



hence E x ((M^)* n ) p < C^nE\fj 1 \ p . (5.8) 
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Now we can conclude similarly to Locherbach, Loukianova and Loukianov 
(2011): For 1 <p < 2, 



JP x (\N t /t-£\>£e) 

< (2 p - 1 lE x \R 1 - l/£\ p ^ + 2 2p+2 C p lE x \f] 1 \ p £) 1 ' 

< c(p)v(x) ((i + (i/ey- ') + m p £c p ) ' 



e p tP- 1 

This concludes the proof. □ 

Once the deviation inequality for the counting process (N t ) t is proven, we 
obtain on the lines of Locherbach, Loukianova and Loukianov (2011), The- 
orem 3.2, the following general deviation inequality for additive functionals 
of the original Markov process X, built of bounded functions. 

Theorem 5.2. Grant Assumptions 2.1 and 2.2 with $(t>) = cv a , < a < 1. 
Put p = 1/(1 — a). Let f G L 1 ^). Suppose that \\f\\oo < oo. Let x be any 
initial point and < e < ||/||oo- Then for all t > 1 the following inequality 
holds: 



i r* 

t 



f(X s )ds-n(f) 



> : } < K(£,p,v,X)V(x)r (p - 1) x 

x i in f hp ifi< P <2 r } 



e p 1 1 J Moo 



Here K(£,p, v,X) is a positive constant, different in the two cases, which 
depends on £, p, v and on the process X through the life cycle decomposition, 
but which does not depend on f , t, e. 

Remark 5.3. The above result holds for any starting measure v such that 
JE U (R P ) is finite, so a fortiori for any measure v such that V G L 1 ^). In 
contrary to most of the existing results in the literature (see e.g. Cattiaux and 
Guillin (2008)) we do not need to assume absolute continuity of the initial 
law of the process with respect to the invariant measure /i. 

Proof. First of all, since the law of X starting from a fixed point x is the 
same as the law of Z 1 starting from the initial measure 1P X , we certainly have 
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that 



f(X s )ds-n(f) 



> e 



iR, 



f{Z])ds-nU) 



> e 



Now for p < 2 the rest of the proof is exactly the same as the proof of Theorem 
3.2 in Locherbach, Loukianova and Loukianov (2011). The only difference 
compared to there is that the variables £ n = J R n+1 (f — ji{f)){Zl)ds are no 
longer independent but only 2-dependent. Hence, the same trick as in the 
proof of Theorem 5.1 applies: one has to separate even and odd terms. But 
this does only change the constants in the upper bound. 

For p > 2, we use the Fuk-Nagaev inequality again. We start as in 
the proof of Theorem 3.2 of Locherbach, Loukianova and Loukianov (2011). 
Denote 

N t 



a 



Put f:=f-fJL{f). Then 



t 



-t 



< £5 } , 5 = e 



< 1. 



IP, 



< JPt 



+JP X 



f(Zl)ds - tii(f) 



> te 



f{Z l s )ds 

R-N t + 1 



te . 
> - I + P, 



l]V t + l 



RzDds 



f{Z])ds 



Ri 

p x (Q t c ; 



>|;a 



= A + B + C + D. 



The terms A and C are handled as in Locherbach, Loukianova and Loukianov 
(2011). Term D is controlled thanks to Theorem 5.1. So the main term that 
has to be controlled is the term B, and we have 

\ k=l J 



Put l k = j±-S k , then 4 = Jl k+1 (g - ii(g))(Zl)ds, where g = f/\\f\U 
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te 



\\g\\oo = 1- We write 

/ IW+S)] \ f [*«(!+*)] 

y k=i J \ k=i 

< Pj sup S k >tS/3], 6 = e/\\f\\ 00 <l, 

\k<[U(l+S)] J 

Sk = Yli-i \£k\, an d apply the Fuk-Nagaev inequality of Theorem 7.1 with 
n = t£(l + 5) and A = t5/ 12. This gives the following upper bound 

/ [«(!+*)] \ 

where 

m p = E(\^\ p ) < 2 p E((R 2 - Rtf), a 2 = Varfa) < 4 E((R 2 - R,) 2 . 

Therefore, since 5 = £/||/||oo, there exists a constant K(l,p, v,X) depending 
only on the process and the life cycle decomposition, but not on the function 
/, such that 

/ [«(!+*)] \ 

JP X \ £ U>te/?>\ <K(lwX)\\f\^%-^h-W. 
This finishes the proof. 

□ 

6 Examples 

We close our paper with two examples where the above deviation inequalities 
can be applied. 

6.1 Multi-dimensional diffusions 

Consider the solution of the following stochastic differential equation in M d 

dX t = b(X t )dt + a(X t )dW u 
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where W t is an n— dimensional Brownian motion, n > d, such that b is a 
locally bounded Borel measurable function M d — » M d and a is a bounded 
continuous function M d — > M dxn which is uniformly elliptic: Writing a := 
aa*, we suppose that there exists e > such that 

<a(aOe,e>>e|K|| 2 

for all x G Classical results on lower bounds for transition densities 
of diffusions (see for instance Kusuoka and Stroock (1987)) imply that in 
this case any compact set of M d is petite. We cite the following recurrence 
conditions from Fort and Roberts (2005). Suppose there exist M,f3, r y > 
and / < 2 such that 

sup \\x\\ ~( 2+ ') < x, a(x)x >= j3, sup \\x\\~ l tr(a(x)) = '-f, 

x:||x||>M x:||x||>M 

sup < b(x),x >= —r, for some r > (7 — j3l)/2. 

x:\\x\\>M 

We choose 

and put m = 2 — 1+k, thus 2— m = 1 — k. Let V(x) = \\x\\ m outside a compact 
set. Then sup x .^ >M AV(x) < 00, and standard calculus shows that for all 
\\x\\ > M, 

AV(x) < m f-r + ^[7 + (m - 2)f3] 
Then by our choice of k, f :— r — ^[7 + (m — 2)f3] > 0. Hence for ||x|| > M, 

AV{x) < -$oV(i), 

where 

2-1 

$(x) = mr x a , with a = < 1. 

m 

Hence we get polynomial moments of regeneration times up to the order 
m/(2-l) = 1 + k/(2-Z). 
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6.2 Solutions to SDE's driven by a jump noise 

This chapter is inspired by a recent work of Kulik (2009) on exponential 
ergodicity for solutions to SDE's driven by a jump noise. More precisely, 
consider the solution of the following stochastic differential equation on M d 
driven by a jump noise 



dX t = b(X t )dt+ / c(X s _,u)n(dt,du) + / c{X s _,u)n{dt,du). (6.1) 
J\\u\\<i ^I|m||>i 

Here /i is a Poisson random measure (PRM) on 1R + x lR q , having compen- 
sator fi(dt,du) = dtv(du), and j2(dt,du) = /i(dt,du) — dtvidu) denotes the 
compensated PRM. We follow Kulik (2009) and impose the following condi- 
tions on the coefficients b and c. The drift function b belongs to C 1 (]R d , M d ) 
and satisfies a linear growth condition. The jump rate c(x,u) is one times 
continuously differentiable with respect to x. Moreover, 

\\c(x,u)-c(y,u)\\ < K(l+\\u\\)\\x-y\\, \\c(x,u)\\ < ^(x)||u||, x,y G R d ,u G R q , 

where K is some constant and where ip : lR d — > 1R + satisfies a linear growth 
condition. Finally we impose a moment condition on the Levy measure v. 
For all R > 0, 

/sup (\\c(x, u)\\ + ||V x c(x, u)\\) v(du) < +oo. 
i:||x||<fl 

Then for any fixed x G M d , there exists a unique strong solution X t to (6.1), 
which is a strong Markov process, having cadlag trajectories. 

We quote sufficient conditions implying that compact sets are petite from 
Kulik (2009). For this sake, we have to introduce some notation. Let S q = 
{v G M q : \\v\\ = 1} be the unit sphere in M g . For any w G S q and for any 
q G]0,1[, let V+(w,q) = {y G M q : < y,w > > g\\y\\}, and V(w,g) = {y G 
M q : | < y,w > | > g||y||}. Then Kulik (2009) obtains the following result. 

Proposition 6.1 (Kulik 2009). Suppose that the following assumptions hold. 

1. Cone condition: For every w G S q , there exists g g]0, 1[, such that for 
every 5 > 0, 

u(V(w,q) H {u : ||u|| < 5}) > 0. 
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2. N on- degeneracy condition: There exists a point rr* G M d and a neigh- 
bourhood of x* such that c(x,u) = x( x ) u + S(x,u), for all x e O*, 
and 

||<J(a:*,u)|| + \\VJ{x*,u)\\ = o(\\u\\), 

as \\u\\ — > 0. Moreover, the functions b(.) = b(.) — f^ <1 c(.,u)u(du) 
and x are one times continuously differentiate and satisfy the joint 
non- degeneracy condition 

rank (vb(x*)x{x*) — Vx(x*)6(x*)^ = d. 

3. Support condition: For any R > there exists t such that for all x with 
\\x\\ < R, 

x* e suppP t (x, ■). 

If the above conditions hold, then any compact set is petite. 

Remark 6.2. 1. In the one- dimensional case d = q = 1, the above conditions 
can be stated in a simpler way. For example, condition 1. can then be written 
as follows: For all 5 > 0, v(u : < \\u\\ < 5) > 0. 

2. Simon (2000), Theorem I, gives a sufficient condition for condition 3. 
above to hold, see also Proposition J^.l in Kulik (2009). 

Proof. Theorem 1.3, Proposition 4.3 and Proposition 4.4 of Kulik (2009) 
show that under the above conditions, the following Dobrushin condition 
holds: For all R > 0, there exists t* = t*(R) such that 

inf [ [P t .(x,-)AP t .(y,-)](dz)>0, (6.2) 

where for any two probability measures P and Q, 

t p A Q ^ - {wTq) {z) a WTq) <*>) (p + Q){dz) ' 

From this the claim follows since (6.2) implies that any compact set is a 
petite set. □ 

It remains to give conditions that are sufficient for the recurrence con- 
dition (2.2), (2.3) respectively. There is a wide range of possible conditions 
and in what follows we restrict attention to a particular sufficient condition 
which is stated in the same spirit as the conditions of Proposition 4. 1 of Kulik 
(2009). 
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Proposition 6.3. Suppose that the conditions of Proposition 6.1 hold. Sup- 
pose moreover that there exist M, 7 > and < I < 1 such that 

1. Moment- condition: There exists m > 1 such that \\u\\ m v(du) < 
00. 

2. Moderate jumps: The function c can be decomposed into c = C\ + c 2 
such that 

(a) ||ci(a;,u)|| < 7||x|| i ||w||, u G R q , \\x\\ > M. 

(b) \\x + C2(x, u)\\ < \\x\\, \\u\\ > 1, ||x|| > M, and C2(-,u) = 0, if 
\\u\\ < I. 

3. Drift- condition: 8wp x -.\\x\\>M ||x|| _<1+ ^ < b(x),x >= —r, for some con- 
stant r satisfying r > 2j \\u\\ m u(du). 

Then there exists M > M such that (2.3) holds with B = {x : ||x|| < 
M }, B petite, V(x) = \\x\\ m and $(x) = cx^ a , where a = ^ < 1. 

Proof. We use the drift condition for the generator defined for all functions 
F in the extended domain of the generator 

AF(x) =< b(x), VF(x) > + / (F(x + c(x, u)) - F(x) - l {M <i} 

J Mi 

< VF(x),c(x,u) >)v(du). 
Applying this to V(x) = \\x\\ m yields for all > M, 

AV(x) =m<b(x),x> \\x\\ m - 2 + [ (\\x + c(x,u)\\ m -\\x\\ m )is(du) 

J\\u\\>l 

+ / (\\x + c(x,u) \\ m - \\x\\ m - m < x, c(x,u) > \\x\\ m ~ 2 ) v{du) 

J\\u\\<l 

< -m ■ r\\x\\ m - 1+l + [ (\\x + c{x,u)\\ m - \\x\\ m )v{du) 

J\\u\\>l 

+ / (\\x + c(x,u)\\ m -\\x\\ m -m<x,c(x,u) > \\x\\ m ~ 2 ) u(du). (6.3) 

J\\u\\<l 
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We start with the term in the last line. By Taylor's formula, writing h = 
c(x,u) = Ci(x,u), since ||w|| < 1, we certainly have that 



I m— 2 



\x + c(x, u)\\ m — \\x\\ m — m < x, c(x, u) > \\x 
< - sup I < h,V 2 V{y)h > \ 

2 y£]x,x+h[ 

<I m [l + |m-2|]||/i|| 2 sup ||y|r- 2 . 

y£]x,x+h[ 

Here, ]x, x+h[ denotes the d— dimensional interval ]x±, Xi+hi[x . . . x]x d , x d + 
h d [. 

Applying condition 2. (a) to h — Ci(x,u), where < 1, yields 

||/i|| 2 <7 2 |N| 2i ||«|| 2 . 

If m - 2 > 0, we choose M > M such that (1 + ^M 1 ^ 1 )™- 1 < 2 (recall that 
I < 1). Then we obtain 

sup |M| m - 2 = \\x + h\\ m ~ 2 < \\x\\ m ~ 2 [l + 7 ||x||'- 1 ] m_2 

ye]x,x+h[ 

< || a ;ir- 2 [i+7M^- i ] m - 2 

< 2||a;|r- 2 . 

If m < 2, we can proceed similarly, 

sup IMP" 2 < \\x\\ m - 2 [1 - ^WxW 1 - 1 }" 1 - 2 

y&]x,x+h[ 

< 2||xir~ 2 , 

where we choose M such that (1 - 'jM l Q ~ 1 ) m ~ 2 < 2. 

As a consequence, for any \\x\\ > M , the last line of (6.3) is bounded 
from above by 

m ([I + \m- 2|] 7 2 / \\u\\ 2 is(du)) ||^||— 2h- 2 z < Q ^-i-i^um-i^ (g 4) 

since < M 1 ^ 1 . Here, M 1 ^ 1 — > as M — > oo, and C is some constant. 

Hence the last term of (6.3) will be neglectable for our purposes. 
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Concerning the first jump term in (6.3) we proceed as Kulik (2009), proof 
of Proposition 4.1: For > 1, using condition 2. (b), we have 

\\x + c(x,u)\\ m - \\x\\ m < \\x + c(x,u)\\ m - \\x + c 2 (x,u)\\ m 

= \\x(u) + c 1 (x,u)\\ m -\\x(u)\\ m , 

where x{u) = x + c 2 (x,u), and then, applying Taylor's formula, 

\\x(u) + ci(x,u)\\ m - \\x(u)\\ m < m||ci(x,u)|| sup IMP -1 - 

ye]x(u),x(u)+ci(x,u)[ 

Now, since m > 1, we argue as before and obtain, using successively condition 
2. (a) and 2. (b) and ||u|| > 1, 

m||ci(x, u) || sup IMP" 1 < m 7|l a; ll'll ti ll (ll x (' u )ll + tIMI'IMI)™ 

ye]x(u),x(u)+ci(x,u)[ 

, m—1 



< m 7 ||a;|| z ||-u|| (||x|| + 7 ||a;||'||-u| 

< m 7 ||a;|r- 1+ '||M|r(l + 7M^ 1 ) m ~ 1 

< 2m 7 ||x|r- 1+ '|| M |r t , 

by the choice of M . As a consequence, the first jump term in (6.3) can be 
upper bounded as follows: 

/ {\\x + c(x,u)\\ m - \\x\\ m ) v{du) <m\\x\\ m ~ l+l 2 7 f \\u\\ m v{du) . 
J\\u\\>l L J\\u\\>l 

Collecting all the above results, we finally obtain that for all ||x|| > M , 
AV(x) < m ( -r + 2 7 ! \\u\\ m v{du) + CM^ 1 ) ^^r r 

By condition 3., for M sufficiently large, — r+2 7 J^ >x \\u\\ m v(du)+C M 1 ^ 1 < 
eventually, and this implies the assertion. □ 



7 Appendix 

For the convenience of the reader we give in this section a Fuk-Nagaev in- 
equality for sums of two-dependent identically distributed centred random 
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variables admitting a moment of order p. This inequality is the key tool to 
our deviation inequalities, and its proof can be found in the book of Rio 
(2000). 

Let Xi,X 2 , . . . be centred identically distributed random variables which 
are two-dependent and such that E(\Xi\ p ) < oo for some p > 2. Put S k : = 
Xi + . . . + X k . Then the following deviation inequality holds. 

Theorem 7.1. For any A > we have that 

P (sup \S k \ > 4A j < C(p) {a 2ip -^\- 2ip - l ^n p - 1 + m p n X' p ) , 

\k<n J 

where a 2 = Var(X 1 ) and m p = E(\Xi\ p ). 

Proof. We use the Fuk-Nagaev inequality presented in Rio (2000), Theoreme 
6.2. First of all, since the variables are two-dependent, we certainly have the 
upper bound on the a— mixing coefficients a n = sup fe>ra a(a(Xi), o~(Xk+i)) : 

a = -, a\ < -, a n = for all n > 2. 

Hence, using the notation (1.21) of Rio (2000), we can upper bound 

a-\u) < 2 l [0i i[( w ). 

As a consequence the expression R(u) of Theoreme 6.2 of Rio (2000) is given 

as 

R{u) < 2Q(u)l [0) i [ (u) < 2Q{u), 

where Q(u) = inf{x : H Xl (x) < u} is the quantile of \Xi\, H Xl (t) = P(\Xi\ > 
t). But since X\ admits a p— th moment, we certainly have that 

H Xl (t) < m p t~ p , 

by Markov's inequality (recall that m p = E(\Xi\ p )). Since Q is the gener- 
alised inverse function of H Xl , this implies that 

Q{u) < m l J p u~ x l p , 

and this in turn leads to 

H(u) = R-\u) < 2 p m p u~ p . 
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Now we can apply (6.5) of Rio (2000). First of all notice that by the two- 
dependency structure 

n n 
i=l j=l 

Thus we obtain, for any r > 1, 

P ( sup 1^1 > 4A ) < 4(1 + — ) + 4nA _1 y Q(u)du 



k<n J \ rs 2 



/ \ 2 \ -r/2 pH(\/r) 

< 4 M + — J + AnX- 1 / ml /p u- 1/p du 

( X 2 Y r/2 4nX' 1 m 1 p /p /w u i 

/ \2 \ -r/2 

< 4 1 + - J +(7(p)m p nA- p r^ 

/ \ 2 \ -r/2 

< 4 ( 1 + ^^j +C , (p)m p nA- p r p ' 1 . 

Here, C(p) is a constant depending only on p. Now we choose r = 2(p — 1). 
By assumption on p, r > 1. Finally we get 

P (sup |S fc | > 4A ] < C{p) (a 2{p - l) X- 2{p - l) n p - 1 +m p n X~ p ) . 

\k<n J 

□ 
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